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1 Introduction 


C. G. J. Jacobi called December 23, 1751 the birthday of the theory of elliptic functions [T], 
but birth is preceded by conception. In this case, conception occurred 33 years earlier, in 
1718, when Giulio Garlo Fagnano published three short papers concerned with measuring 
the length of an arc of the lemniscate [2]. Among other things, Fagnano discovered how to 
double the arc: 

“Given the two equations 

= Vi-v^ (1) 

vl - Z 

dz 2du . 

I say that if the first holds, so does the second.” 

Carl Ludwig Siegel adopted Fagnano’s work as the starting point for his three volume 
introduction to topics in complex analysis [ 6 ] . He proposed a reconstruction of how Fagnano 
might have been led to his proof, and also simplified the argument, making it more elegant, 
and transforming the formulae into a more modern form, which is 


r 


2 


2iv? 

1 — ’ 



(3) 


Here we recognize an example of the complex multiplication of elliptic curves. David Hilbert 
said that the theory of complex multiplication of elliptic curves is not only the most beau¬ 
tiful part of mathematics, but of all science0No wonder Siegel chose to begin with it. 
The substitutions r = \/^, u = \/^ convert eq(|^ to 

2 s /-^ dt /-^ ds 

^ ~ / a 2 1 ’ / / , ,o —(! + *) / / . o (4) 

V Js^ — 1 Jo V — t Jo V 4s^ — s 


This leads us to the WeierstraB p-function as the uniformizing variable, which satisfies the 
differential equation p'^ = 4p^ — 92 P — do p] here, g 2 = 1, 53 = 0. Hence the Klein modular 
invariant j(T) has the value j{l J- i) = j{i) = 1 , so the lattice of periodicity D for p{z) is 
generated by 1 and i. Inverting the integrals, we find an expression for p ((1 J- i)z); this is 
a complex multiplication of the Weierstrass p-function. Amongst other things, it implies 
that (1 J- z)H C D. Equation Q yields the complex multiplication formulsj^ 


p{il + i)z) 


8 V P(^) J 


1 

2 

3 


Fagnano’s words, rendered by Ayoub [3- 

In 1932, at the International Congress of Mathematicians. [3], p.2 00. 

We shall recover this formula in a different way below. Cp. eq(46l et seq. 


( 5 ) 
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Every elliptic function is doubly periodic with respect to a lattice in C of rank 2, and 
every such lattice 12 admits multiplication by a rational integer. That is, for n € Z*, nO. C 12. 
Lattice-invariant functions express this property through a real multiplication formula. For 
instance, for any lattice, the duplication formula for p{z) corresponding to the inclusion 
212 C 12 is 

which, after p" and p' are expressed in terms of p by means of the differential equation, 
yields a formula that expresses p{2z) as a rational function of p{z): 


P{2z) 


1 (12p(z)^-ff2)^ 

16p(2:) (4p(2;)2 - 52 ) 


2p{z). 


( 6 ) 


This situation is reminiscent of the duplication formula for sinz, which can be expressed 
in four ways, viz. 


sin 2z = 


Although every non-zero rational integer is a multiplier, it is unusual for a lattice to admit 
a complex multiplier, a complex number /r such that //12 C 12. The necessary and sufficient 
condition is that 12 be the ring of integers of an imaginary quadratic algebraic number field, 
and then /r can be any non-zero integer in the field. In the lemniscatic case, 12 is the ring 
of Gaussian integers Q[i] in the field Q(i). 


2 sin z cos 2 ; 

(2sin(2;) 


2 sin 2 


dz 


2 sin 2 \/1 — sin^ 2 ; 

2 sin 2 ; sin(z -|- 7r/2). 


*** 

Superficially, multiplication formulae seem to arise as limiting special cases of an ad¬ 
dition formula, and the addition formula appears to emerge from the periodicity of the 
functions with respect to an underlying lattice. From a more general perspective, the ex¬ 
istence of a multiplication formula is a consequence of a theorem of Siegel that asserts 
the algebraic dependence of n -|- 1 meromorphic functions on an n-dimensional compact 
complex manifold [5]. Here, the manifold is the elliptic curve, that is, the complex torus, 
on which the meromorphic functions live, and C is its universal covering space. Given the 
algebraic dependence theorem, it follows that f{z) and /{/iz) satisfy an algebraic equation 
for any meromorphic function / defined on the manifold and any complex number ^ that 
carries the lattice of periods induced on the covering space into itself. 
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There is no a priori reason, however, to believe that multiplication formulae do not exist 
in situations where meromorphic functions do not live on a compact complex manifold. Of 
course, some structure must be presupposed in order to prove anything. We have in mind 
inflationary tessellations - particularly, aperiodic tessellations of the kind exemplified by 
Penrose tilings - as the geometric substratum that substitutes for the lattice of periods in 
the classical theory, and functions defined in a way that mimics definitions of some of the 
most important classical higher transcendentals. 

One purpose of this paper is to construct multiplication formulae using the properties 
of positional representation for numbers. When the radiji[^ is a positive integer, the results 
recapitulate multiplication formulae for many classical functions. If the radix is not a 
rational integer, and particularly if it is a complex number that is not a quadratic imaginary 
integer, both the functions and the multiplication formulae appear to be new. 

Although these results may be interesting, they cannot be said of themselves to be 
important. Some of the new functions are similar to functions that appear in the the¬ 
ory of elliptic functions and higher arithmetic, but analogues of the properties that make 
those functions significant are thus far unknown. What may make the new functions, and 
their multiplication formulae, of some immediate interest is their connection to inflationary 
tilings and particularly to those that are aperiodic. 

This observation sets the agenda for the paper. We make use of the connection between 
inflationary tessellations and positional representation described in a previous paper [3] to 
connect these tilings to positional representations, which can then be used to construct 
meromorphic functions on C that incorporate properties of the tiling structure. If, for 
instance, the inflationary tiling is a non-degenerate lattice in C (hence, periodic), then the 
functions are elliptic functions. 

Tessellations built up from an unlimited number of types of tiles are uninteresting. As 
was shown in [3], if an inflationary tiling employs a finite number of types of tiles (the 
Penrose tiling uses 2), then the inflation factor - which turns out to be the radix of the 
associated positional representation - must be an algebraic integer. Thus, although multi¬ 
plication formulae exist for any radix, it is only the ones that are algebraic integers that 
couple function theory to an inflationary tessellation. This is another connection with, and 
an analog of, complex multiplication. We shall begin by deriving various unrestricted gen¬ 
eralizations of classical multiplication formulae, but the heart of the paper concentrates on 
radices that are algebraic integers related to the golden number (j) = In particular, we 

explore the radix It satisfies + x^ = 1 and is the multiplier for a particularly simple 
aperiodic tessellation of C that uses four types of tiles, three of which are geometrically 
similar. We call this the special tessellation. 

^ Often called the ‘base’ of the number system. 
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The paper is organized to develop in detail this specific example of an aperiodic tiling 
of the plane and a family of associated functions that have complex multiplication formu¬ 
lae of the new kind. But first we examine some classical functions from the standpoint 
of positional representations for numbers in M and in C. Section introduces positional 
representations, both real and complex, in a form appropriate for extension to inflationary 
tilings. Section introduces the positional representation approach for some higher tran¬ 
scendental functions, and applies it to derive generalizations of some classical dnplication 
formnlae. When the radix is 2, these generalizations specialize to the classical dnplica¬ 
tion formulae. In particular, these ideas are applied to 1-dimensional real mnltiplication 
formnlae such as eq(|36[), for example, which generalizes the classical duplication formula 
for TTcotvrz, a trigonometric function that is strncturally a 1-dimensional version of the 
Weierstrafi elliptic function p(z). 

Positional representations fall into two main classes: those for which the ‘integers’ of 
the representation relative to the radix p have a unique representation, and those which do 
not. In section]^ we investigate two examples for which representations of p-integers are 
not unique. 

One of the pleasant peculiarities of our approach is that partition of integers according 
to the number of digits in a positional representation leads to exact formulae for some 
associated functions. In section exact identities for partial sums of certain meromorphic 
functions defined by infinite series are derived. Since the partial sums are finite, questions of 
convergence are (temporarily) avoided. The Riemann (^-function has a ready generalization 
to this setting. Section considers the generalization when the radix is the golden number 
(j) = and develops a formula that provides an analytic continuation to the complex 

plane. 

Next we make explicit the connection with aperiodic tilings. Sectionj^applies these ideas 
to a ‘duplication’ formnla for an aperiodic (1-dimensional) tessellation of K associated with 
the radix —(j). This is preliminary to the work of section which takes up the analogous 
problem for the special tessellation, an aperiodic tiling of the plane associated with radix 

Section examines the convergence of the previously introduced infinite series. The 
paper concludes with remarks collected in section [10} 


The principal specific resnlts of the paper are for the complex radix p = presented 
in section 8.2 the exact duplication formula for partial sums of the p-Hurwitz ((-function 
in eq(92), and the corresponding identity for the Weierstrafi p-function pp{z) in eq(94). 
The poles of Pp{z) are the vertices of the special tessellation. 

Another formnla of interest provides an analytic continuation for the p-Riemann ^- 
function, which specializes to the usual Riemann ^-function and appears to provide some¬ 
thing new even in this case; cp. eq(l^. 
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2 Positional representations 


Many classical examples of multiplication formulae can be derived from the following 
paradigm. Given a fixed integer p > 1, express an arbitrary non-negative integer n in 
the form kp + r with r = n mod p. Terms that have the same remainder mod p in a sum or 
product representing a function are collected. This, apart from questions of convergence, 
supplies the desired formula. 

The positional representation of n in radix p also expresses these relationships. If the 
digits of the representation are the integers A = {0 ,l,...,p—1 }, then 

n = 4 G A, 

k>0 

where the sum is finite. This procedure can generalized. The radix need not be positive, nor 
an integer. Any p G C such that |p| > 1 can serve as a radix for a positional representation. 
The set of digits can be selected in many ways, but the number of digits is [ ], 

where dim(p) = 1 for representation of real numbers and dim(p) = 2 for representation of 
complex numbers [x] is the least integer greater than or equal to x. 

In either event, for z running through some set of positive measure (which may be all 
of K or C), 

z = ^ZkP~’^, Zk^l^, (7) 

where finitely many positive powers of p appear in the sum. The series converges absolutely 
since it is majorized by a geometric series. 

It will be technically convenient, and easier to express the ideas we want to emphasize, 
if we limit ourselves ourselves to radices for which the number of digits is 2. In this case, 
A has two elements, which we select to be 0 and 1. Then the digits Zk in eqQ are bits. 
This case already captures the essentials, both of complications and of consequences, of an 
arbitrary radix. Since only two digits occur, the corresponding generalized multiplication 
formula will be referred to as a p-duplication formula. 

The set 

72:= , 4gA, (8) 

is naturally called the set of p-remainders for the positional representation. It is R, and 
certain natural subsets of R, that will appear as the representative tiles of the associated 
tessellation. 

Let A[A], resp. A[p], denote the set of polynomials in X, resp. p, with coefficients in 
A. The elements of A[p] are called p-integers because A[p] is the usual set of non-negative 
integers when p is a positive integer greater than 1. 


®Roughly speaking, dim(p) is the dimension of the remainder set R defined by eqj^ below. 



Let A„[X], resp. A„[/9], be the subsets of polynomials in X, resp. polynomials evaluated 
at /?, that have n coefficients, the n-digit polynomials. The leading coefficient may be 0. 
The multiplicative containment XA„[A] C A„_|_i[A] is the essential ingredient for deriving 
multiplication formulae. 

An, [A] has 2"' distinct elements but the number of elements of A^fp] that have distinct 
values may be less. The coincidences are the consequence of algebraic identities satisfied by 
p. Suppose that p{X) and q{X) are relatively prime polynomials with coefficients from A. 
If p{p) = q{p), then p and q are distinct representations for the same number. In particular, 
we see that in order for coincidences to occur, p is necessarily an algebraic integer. 

Let us write p ~ g if p{p) = q{p). is an equivalence relation. The collection of 
equivalence classes A[A]/(~) is in 1-1 correspondence with the values assumed by the 
polynomials. 

Specific sets 12 of representatives of A[A]/(~) play an important role in what follows. 
The important point is that a set of unique representatives - not necessarily complete - of 
p-integers be available. 

It is no accident that we use a notation - 12 - that was previously used to denote a 
non-degenerate lattice in C, because when the lattice has a complex multiplication, it is a 
set of unique representatives of integers of the imaginary quadratic number held associated 
with the complex multiplication. 

It was shown in [1] that the inhation factor - the “multiplier” ~ for an inhationary 
tiling must be an algebraic integer, whence it follows that in the applications of positional 
representation to tessellations, some p-integers may have more than one positional rep¬ 
resentation; in other words, the equivalence relation will not be trivial]^ In the realm of 
p-multiplication formulae for meromorphic functions, this has the consequence that some 
terms of inhnite series summed over the set A[p] will be repeated sufficiently often to force 
divergence. The obvious solution is to select a unique representative polynomial from each 
equivalence class, i.e., to sum over the elements of 12. This selection will depend on the 
properties of the radix. 

Given 12, let 12* denote the set of non-zero elements, and denote the subset of n-digit 
p-integers by f2„. Elements of I2„ may have 0 as the leading digit. It will be helpful to have 
a notation for the n-digit p-integers whose leading digit is 1, namely 

Xn •— I2fj n > 1. (9) 

The set I2i of 1-digit p-integers is exceptional; we put I2i = Vi = {0,1}. The V^, are 
pairwise disjoint, and 12 = U^i Xn- Put 12* = 12^ — {0}. 

®We cannot claim that there will necessarily be coincidences because the minimal polynomial of p may 
not have coefficients from ±A. Thus the representations of p-integers for p = y/2 are unique. 
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3 The general paradigm 


As an introduction to the method, we shall formulate the general approach and then 
specialize it to obtain examples for classical functions. Afterwards we will head into the 
uncharted waters and unknown functions associated with radices for which the represen¬ 
tations of p-integers are not unique. Throughout, we assume A = {0,1}. 

3.1 The p-Hurwitz (^-function 

Assume that = 2. Suppose that p{p) / q{p) for distinct polynomials p,q € A[A], 

so that n = A [p]. 

The condition of uniqueness - i.e., that each p-integer is represented exactly once in 
A[p] - is satisfied if p is a positive integer. It is also satisfied if p is transcendental, for then 
the only solution to p{p) = q{p) \s p = q. 

Algebraic cases sometimes lead to uniqueness. Here are two examples: 

(1) If p = 3/2, then the only solution to p(p) = q{p) \s p = q. Suppose otherwise. After 
eliminating common powers of p and clearing the denominator, the leading term is a power 
of 3, and the trailing term a power of 2. 

(2) If p = y/2, then a p-integer can be written in the form m + ny/2 with m and n 
non-negative integers. Again, p = q is the only solution. 

Cases where the representation of integers is not unique are of greater interest because 
they may be related to tessellations. For this reason we shall explore p = (j) = , the 

‘golden number’, and p = i\/(p. Here (j)'^ = cj) + 1. The ratio of unique n-digit (/-integers to 
the number of polynomials in A„[A] approaches zero as n increases. Nevertheless, we can 
specify a subset of polynomials such that every (/-integer has a unique expression, and let 
the sum in an infinite series run over that set. This will be done in section HI 

With Q, denoting a set of unique representatives of A[p]/(~), let denote the subset of 
H of representatives of A„[p]. Set Q* = Cl — {0} and — {0}. Where the p-positional 

representation of integers is not unique, we will explicitly state how the elements of Cl are 
chosen. We will write #{S) for the number of elements in a set S. 


Define the p-Hurwitz (/-function. 


Cpis,z) ■= 


1 

(z -b ojy ’ 


and the p-Riemann /^-function. 


Cp(s) := 


1 

UJ^' 


( 10 ) 


( 11 ) 
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Discussion of convergence of the series will be postponed to section because it depends 
on particular properties of the radix. 

Suppose that the series converge absolutely and uniformly on compact subsets of the 
complement of D, resp. D*, and that the representation of p-integers in A[p] is unique. 
Then 


p"Cp{s,pz) 


{z + uj/pY 


E 

y_^ 

(z + w)' 

t<;er2 ' a;en 


+ E 


{z+l+wY 


Hence 


P^ Cp{s, pz) = Cp{s, z) + Cp(s, z + 1/p), 


( 12 ) 


which is the p-duplication formula for the p-Hurwitz (^-function. 

The function Cp('S, ^ + l/p) is holomorphic in a neighborhood of z 
find 


0. From eq(12) we 


Cp{s,l/p) = Y^^{p" C,p{s,pz) - Cp{s,z)) 

1 \ 1 

^ {pz+ 

= {p^-ycp{s). 

Hence 

{p^-l)Cp{s) = Cp{sMp) ( 13 ) 

expresses the p-Riemann (/-function in terms of the p-Hurwitz (/-function. We have implic¬ 
itly used uniqueness to conclude that — 1/p is not a pole, for if it were, then it would have 
the same value as some cu G D* and clearing the denominator would yield two representa¬ 
tions for 1. 

The poles of Cp(s, z) are the points — flj^We shall calculate the Laurent series for Cp(s, z) 
at z = 0. In a neighborhood of z = 0, the function 

h{z) := Cp{y -= 7-^ 

^ ^ ^ {z + ujY 

^It would be aesthetically more satisfactory to define the Hurwitz function, etc. so that the set of poles 
is Q, but, for the moment, we prefer to make the connection with classical notations. 
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is holomorphic. Then 


Therefore 


d^hjdz^ 


^ ( -iNfc r(s + fc) ^ 1 

-0 ^ ^ r(s) + 


2 = 0 




(p(^, 2) — —+h-i .h C(i(^+ 

^ k=0 


r(s)r(A: +1) 


(14) 


Bearing mind that the constant term on the right side of eq(14) is Cp{s), evaluate both 
sides at z = 1/p and use eq([l^ to find 


Cpis) = 


-2 


i+E(-o 


fc>i 


k r(s + k) Cp{s + k ) 
T{s)T{k + 1) 


(15) 


Equation © provides an analytic continuation of the p-Riemann ^-function to the entire 
complex plane. 

If eq(|15[) is valid (i.e., for values where the series converges appropriately), the function 
has a simple pole at 

s = log 2/log p. (16) 

3.2 Classical examples: p G K, A = {0,1} 

Classical examples of duplication formulae fall into several classes. For a real radix, |p| = 
n G Z"*“: that is, an endomorphism of the lattice Z. The ancillary condition [|p|] = 2 implies 
p = ±2. For p G C, the radix is an endomorphism of a non degenerate lattice in C and 
hence a rational or algebraic integer in a quadratic imaginary number field that depends on 
the lattice. In this case, the ancillary condition is |pp = 2. There are just three essentially 
different cases, which we shall consider separately: p = i\/2, p = , and p = (—1 + i) 

We begin with p = 2. 


3.2.1 The Riemann ^-function 

When p = 2, the p-Hurwitz and the p-Riemann functions specialize to the classical Hurwitz 
and Riemann (^-functions, respectively. 

Equation ( |T^ yields Cp{s) = Cp(S) 1) = C{s), where (/{s) is the Riemann ("-function, 
which expresses YlT=o (2fc+i)‘» ™ terms of ((s). From eq(16) we find the pole of ((s) at 
s = 1. 


®These multipliers are only determined up to a unit of the imaginary quadratic field. 
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( 17 ) 


Equation (|15[) reduces to an identity for the Riemann ^-function, namely 


C(s) 


A I .Nfc r(s + A;) C(s + k) \ 
r(s)r(A: + l) 2^+*^ j' 


which is probably known but I have not found it in the literature. If s = ^ + R is a zero, 
it implies the curious property 

E C (k + 1/2 + it) = - 2 l/^« r(l/2 + ,t). ( 18 ) 

Moreover, if we temporarily write 


\k r(s + fc) ({s + k) 


then, since the zeros of 2® — 2 are 1 + n G Z, it follows that g (^ + 


(19) 


= 0 . 




It may not be amiss to devote a few words to what eq(|17[) means in practice. The series 
converges for Re (s) > 1 and the ("-functions that appear in it are {((s +A:) : k G Z+}, so we 
can calculate ((s) on the critical line Re (s) = 1/2 from the series. Figure[^shows the graph 
of |((^ + iv)\ calculated two ways: the thick pink curve was computed by Mathematica® 
using its internal code; the black curve was calculated using 80 terms of the series from 
eq@. 


Figure 1: Graph of the absolute value of the Riemann (/ function on the critical line using 
analytic continuation provided by eq@ - the black curve - compared with the standard 
function. 

ifi 
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3.2.2 The /9-polygamma function 

The gamma function of Euler has the infinite product expansion 


1 

fM 


OO 



( 20 ) 


where 7 is Euler’s constant. Gauss’ multiplication formula, 

n—1 

r{nz) = ( 27 r)^n^"’-i E 

k=0 


k 

z — 
n 


n 


G Z+, 


( 21 ) 


follows by splitting the infinite product into factors according to the remainders modn. 
The duplication formula is the special case 

r(2z) = ^r(z)r(z + 0 . (22) 

We shall retrieve this duplication formula in the next subsection. 

The polygamma functions, that is, the derivatives of logr(2:), have the expansion 

dz 

(„) _ J-+Mogr(.) 

^ &”+■ 

Evidently {z) is proportional to the 2-Hurwitz (^-function, 

= {—)'^^'^n\C,2{n + l,z), n G (25) 


-7-E 

2: 

fc =0 

OO 


1 


1 


z + k k 

1 

{zTW^ 


(23) 

, 0 < n G Z+. (24) 


so it too will inherit a p-duplication formula. 

Introduce the /9-polygamma function V’p’^E) n G by 

V-W(2) = (-l)-+in!Cp(n + l,z). 

Regarding we must reason somewhat differently. Introduce 


Cp(l,^) 



(26) 


(27) 
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Formally, 


pC,p{'^,pz) = 

= 


Z —' \Z + Lu/p {jjjp 

cjGfi* ^ ^ ^ 


1 1 


+ E 


z ^ vz + w U)J ' ^^\z + l/p + uj l/p + uj 

loGU* ^ ^ uiGfl ^ 


. / X 1 1 

Cp(l) H- -zr-, - — —I" 


E 

ojefi* 


z + l/p 1/p 

1 


1 1 

- + - 


1 


Z + l/p + OJ OJ Ul l/p + u 


Z + l/p^ ^\z + l/p + UJ u 


E 


1 1 


1 


^/p jfe‘Vl/p + ^ ^ 


PCp(l,/5^) = Cp(l,^) + Cp(l,^ + l/p) -Cp(l,l//5), 


(28) 


which is the ^-duplication formula for ^^(1,2;), and hence is proportional to that for 'ij;f'\z). 
The duplication formula for 'ip^f^\z), n > 0 can also be derived by differentiating this 
formula. 

Specializing to p = 2, we find C2(l/2) = 2 - Ylt=i fc(2fc+i) = ^ log 2 so we recover the 
duplication formula for the classical polygamma function: 

2V’^°^ {2z) = 2 log 2 [z] + 'i/^^\z + l /2). (29) 


3.2.3 The p-Gamma function 

Working backward, integration followed by exponentiation of eq@ leads to an infinite 
product we shall call the p-Gamma function: 


1 

W 




n 


1 + 




u 


(30) 


The constant of integration 7p generalizes Euler’s constant, to which it specializes for p = 2, 
in which case V2{z) = r(z). 

From the p-duplication formula eq(28) we find the p-duplication formula for Fp(2;): 


Fp(pz) = e-^^°’(V^)"Fp(z)Fp(z + 1/p). 


(31) 
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3.2.4 The p-cotangent function: p = —2 

The Laurent series for the cotangent function and its derivatives are 

1 1 ' 




TT cot TTZ = - , / , , , 

Z \z + k k J 

fceZ-{o} 


dJ^TT cot vrz 


Etttt 


fceZ 


{z + k)^+^' 


m > 0. 


(32) 

(33) 


The derivatives of tt cot nz can be expressed in terms of the p-Hurwitz ("-function with 
p = —2. Note that here the summation index set Q. of (—2)-integers is Z. To see this, write 


^Ski-2)\ 6k € A 


k>0 


k>0 k>0 

{hk — 2 62 k+i) 4 ^ 


fc >0 


efc4^ : Cfc G {—2, —1, 0,1} 


fc >0 


= Z, 

because {—2,—1,0,1} is a maximal set of integers inequivalent modulo 4 and therefore 


suitable as a set of digits for a radix 4 positional representation. Comparison of eq(33) 
with eq([Io|) shows that 


(-)^d"^7rCOt TTZ 

re! dz'^ 


= (_2(m-I-1, z), m G Z"' 


(34) 




The duplication formula for vrcotTr^; is obtained in the usual way, by integrating the 
series for its first derivative and selecting the constant of integration to force convergence. 
This also works for the p-cotangent function cotp z, which we define as follows, omitting 
the factors of vr for notational simplicity]^ 

Suppose the series 

cotp(2;) ^(35) 


^ \z + (jj CO 


^ThuS COt-2{z) = TTCOtTT^. 
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converges absolutely and uniformly on the complement of for the given p € R with 
1 < |/o| < 2. Recalling that 

Lj = 6 + pQ, 6 G A] uj,uj G A[/5] = n, 


calculate 


p coi p{pz) = - + X] 




1 


1 






coipiz) + ^ 


z + ujIp w/p 

— - + E 

Z + UJ UJ J ^ 
1 




1 


1 




z + l/p + w l/p + ca 

1 


= cotp(z) + 
1 

Wp~ 


z + l/p 

-E 


Z + l/p + UJ l/p + UJ 

+ E' ' 

ojSO* 

1 


1 

Z + 1/p + OJ OJ 

1 


1/p + a; OJ ^ 

= cotp( 2 ;) + cotp( 2 ; + 1/p) - cotp(l/p). 

Thus, with conditions that insure convergence, the p-duplication formula for cotp( 2 :) is 
pcotp(p 2 ;) = cotp(z) + cotp( 2 : + 1/p) - cotp(l/p). (36) 

For p = —2, cot_ 2 (-z) = TTCotTTz and eq(p^ is the classical duplication formula. 


3.3 Classical examples: p G C — K, A = {0,1} 

Now we turn to complex radices. It has already been noted that there are three lattices that 
admit complex multiplications p for which |pp = 2. The associated quadratic imaginary 
field is Q(p). The lattice generators can be taken as {l,p} with p G {i\/2, + *}■ 

The radix p = — 1 + i corresponds to the lemniscatic complex multiplication discussed in 
the introduction!^ In each case the modular invariant is a rational number: 

j(^^/2) = (5/3)3, ■ =-(5/4)3, + + (37) 

^°The unit —1 G Z makes a difference in positional representations with real radix. For instance, A[2] = 
U { 0 } whereas A[—2] = Z. An analogons statement is true for complex radices. The units for 1 + i are 
{±l,±i}. Radix 1 + i yields a snbset of the Gaussian integers Q[z] which, together with i — Q[i], exhausts 
all the Gaussian integers, whereas radix —1 + i produces every Gaussian integer: A[— 1 + i] = Q[i]. 
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3.3.1 The p-Weierstrafi p-function 

Recall the definition of the p-Hurwitz (^-function eq@, 

and the p-duplication formula eq®, 

pXp{s,pz) = Cp{s,z) + Cp{s,z + l/p), 
which is valid for p G C when the series converges. Recall eq@, 

Cp{s,'i-/p) = - i)Cp(s). 


Introduce the p-Weierstrafi function 


1 


Pp{z) ;= + ^ 


1 


1 


Wz + ti;)^ (a;)2 

oj&n* ^ ^ ^ 


(38) 


1 


and its derivatives, 

p<‘)(.) := + 

= {k + 1)! Cp(^ “1“ 2, z), A; G Z^. (39) 

Assuming there are no duplications in A[p], repetition of the argument that led to 


eq(36) for the p-cotangent function yields the p-duplication formula for pp{z): 

P^ Ppipz) = Ppiz) + Pp{z -h 1/p) - Pp(l/p). (40) 

Regarding the Laurent expansion of Pp{z) — l/z^ at z = 0, , the constant term is 0 
(= J2n* (w~^ ~ w“^) ) and eq(39) implies 


d’' {p pjz) - z 2) 
dz^ 


— + !)• Cp{k + 2 ) 


z=0 


SO 


Ppi^) ~ ^ + ^(“)^(^ + l)Cp(^ + 2) z’^. 


(41) 


k=0 


This formula shows that the function (p{k + 2) is the value of an Eisenstein series if p is a 
complex multiplication, and otherwise a generalization of it. 
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When n is a lattice and C O, pp(z) coincides with the classical WeierstraB func¬ 
tion p[z\^). There are no convergence issues. With the restrictions we have placed on p 
(including A = {0,1}) there are, as previously asserted, just three cases. In each case, 
a fundamental domain for 17 C C is the remainder set R = {Y^^=i^k /£ A} and 
the collection {cu -|- i? : w G 17} is a periodic tiling of C by copies of the tile R, and the 
corresponding pp{z) is an even function of z. 

The new cases, of course, will be those for which p is not a quadratic imaginary integer 
of norm 2; the p-integers do not form a lattice; and the associated inflationary tiling is not 
periodic. First let us look at the three quadratic imaginary lattice examples in more detail. 


3.3.2 p = iV2 

For p = i^/2 the remainder set is the rectangle 

R = + iy : Q < X < 1, ^ < y < -y/2| ; 

thus R is the closure of a fundamental domain for the lattice 17 generated by {l,i\/2}. 

This positional representation and R can be expressed in terms of an integer radix 
because z = {iy/2)~^ can be rewritten as 

z = J](<54fc(-2)-2^ + <54fc+2(-2)-(2^+^)) + 

i \/2 ^ + 54 fc + 3 (- 2 )-( 2 '=+^)) 

= x + iV2y, 


where x, ?/ G K are arbitrary numbers expressed in radix —2. The units in Q(i\/2) are ±1. 
The above argument shows that the p-integers are the same for radices ±i\/2- 

The p-duplication formula eq(|40[) specializes to the complex multiplication formula 

- 2 p{iV2z ■. n) = p{z ■. Vt) + p {^z - iV2/2 : 17^ - p (^-7\/2/2 : 17^ . (42) 


The numbers oji = —1 and 0 J 2 = —iy/2 are also lattice generators, and their sum is 
ws = — 1 — i\/2. Choosing z = —1/2, we obtain an identity for the p-function evaluated at 
half periodsp] 


-2p{u2l2) = p{u:i/2) + p(w 3/2) - p(a;2/2). 


that is, 


p(u;i/2) + p(a;2/2) + p(w3/2) = 0. 


This is a standard, and fundamental, identity for the roots of p\z\ 
drop Q. to simplify the notation. 
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3.3.3 

p = ^"*" 2 ^ can be treated similarly. 

3.3.4 p = —1 + i : the lemniscatic case 

The lattice generated by 1 and 1 + i is the ring Q[i] of Gaussian integers; 1 and —1 + i 
also generate Q[i], so the theory of elliptic functions is the same. From the perspective of 
the associated positional representation, there is a difference. A[1 + i] is not the full ring 
Q[i]. One finds A[1 + i] U (i — A[1 + i]) = Q[i] and A[1 + i] 0 (—A[1 + i]) = 0, whereas 
0 = A [—1 + i] = Q[i]. In both cases, the remainder set is the so-called “twindragon” fractal 
region; it tiles A[1 + i] and also C. 


Figure 2: Decomposition of the Gaussian integers Q[i] (in a square neighborhood centered 
on the origin) as the disjoint union A[p] U (i — A)/?]) for radix p = 1 + i. The blue points 
illustrate A[/ 9 ]. 



At this point it will be convenient to make a change of notation. In the defining eq(38) 
for pp{z) we wrote (z+cu)^ where it is customary to write (z —in order to preserve com¬ 
patibility with the classical notation for the Hurwitz zeta function, and the p-generalization 
of it. This notation has the consequence that the poles of pp{z) do not in general coincide 
with the set of p-integers. It made no difference in the previous discussion of p(z) because 
—Q[i] = Q[z]. But now, and for our later work, it makes a difference because we will be 
considering aperiodic tessellations, so it will be preferable to define 


Pp(z) 




In this notation, the p-duplication formula becomes 

Pp{pz) = Pp{z) + pp{z - 1/p) - pp(-l/p). 


(43) 


(44) 
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For p = 1 -\- i the series for pp{z) is a 2-dimensional analogue of the series for the 
polygamma function which runs over the non-negative integers. ^jJ^^\z) is a kind 


of ‘one-sided’ cotangent function (cp. eq(27)). Similarly, pi^i{z) is a ‘one-sided’ Weierstrafi 
p-function. The (1 -|- i)-duplication formula is 

1 -A / 1-i' 


2i pi+i((l -h i)z) = pi+i{z) + pi+i z - 


Pl+i 


(45) 




With p = ( —1 -|- z) we return to the realm of ordinary elliptic function theory and 
complex multiplication. The set of p-integers is the lattice A[—l-|-z] = Q[i], and = 

p{z-.m) = p{^)^ 


Observe that Yjp = ^ is a zero of p[_f_^(z) (it is congruent to a half-period). The 
p-duplication formula reduces to 


- 2i p((-l + i)z) = p(z) + p{z + 


1 + i 


P 


l + i 


(46) 


This is easily transformed into a formula with complex multiplication (1 -|- i) for p{z). 
Substitute z = it to find 


(-2i)p((l -h i)t) = p{it) + p{it + 


1 + i 


P 


l + i 


2 J " V 2 

and use the defining series to show that p{it) = —p{t). Hence eqp^ becomes 

// v\ /\ f 
2i p((l i)z) = p{z) -h p U + —-— J + p 

Now substitute z = 1/2 to find 


(47) 


that is, 


2<>i1^)=(^())+p(5)+p(1/^ 


2 i 63 = ei + 62 + es = 0, 


whence p (^) = 0 and eq(47) assumes its final form. 


2i p((l -h i)z) = p{z) + p{z + 


-l + i 


(48) 


This brings us back to Fagnano’s complex multiplication formula, which can be obtained 


from eq(48) by applying the addition formula for p{z) to the second term on the right, and 
recalling that p'((—1-|-f)/2) = p((—1-|-z)/2) = 0. We immediately find Fagnano’s formula 
eq®: 

p{{l + i)z) = -- ' ' ' ' 


8 V Piz) 


^Again we drop the lattice Q[i] from the notation for simplicity. 
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4 Examples with non-unique representation of /?-integers 


Let p and q be polynomials in the indeterminate X with coefficients in A, and suppose 
that p ^ q. p{p) and q{p) are p-integers. If p{p) = q{p) then the two polynomials represent 
the same p-integer, and this equality is an algebraic equation for p. This shows that non¬ 
uniqueness of the representation can only occur for certain algebraic values of the radix. 
Moreover, every pair of representations of the same p-integer leads to an algebraic equation 
which has p as a root. 

In our case, 1 < |p| < 2 implies that the simplest examples arise as roots of the equations 

+ I, 

so that p = ^(j). Using positional notation relative to the radix, the first case implies the 
identity 100- = Oil-; the second, 110- = 00I-, where denotes the radix point. Each 
equation is the underlying constituent of every non-unique representation of an integer for 
the associated radix. It can be used to construct a unique normal form for each integer. 
This construction will be taken up in the next sections. 


4.1 Dimension 1: Properties of radix (j) = 


Here is a summary of what we know about this radix so far: The golden number (p = 
— 1.61 is the larger root of = x -|- 1. The digits for the positional representation 
with radix (p are A = {0,1}. The (/>-integers are the polynomials in (p with coefficients from 
A, that is, A[(/)]. Let X denote an indeterminate. There are 2” polynomials of degree less 
than n in A [A] but the specialization X ^ (p reduces this number because of the identity 
cp"^ = (p + 1. In positional representation, this identity is 100 = 011 A simple consequence 
is 


L2n 


= i + E' 


>2k-l 


n G Z+. 


(49) 


k=l 


There are two particular expressions that are useful ‘normal forms’: in the maximal 
normal form, the sequence of bits does not contain successive instances of ‘1’; in the minimal 
normal form, the sequence of bits does not contain successive instances of ‘O’. For instance, 
the left side of 1001 = 0111 = 111 is in maximal normal form; the right side, in minimal 
normal form. The number of bits in the maximal normal form is at least as large as in the 
minimal normal form. 

The representation of a (^-integer in either normal form is unique. We shall use maximal 
normal form to construct a list of all (/>-integers without repetition. 

Every binary sequence that expresses an element in A[())] can be written in maximal 
normal form. Suppose 6 = 6n-i ■ ■ - Sq is an arbitrary sequence of n bits with 6 n-i = 1- 
Suppose that 6 kSk-i is the leftmost pair of successive ‘I’s. If A; = n — 1 then 5 is equivalent 


is conventional to drop leading zeros, as well as the radix point when denoting integers. 
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to 1005n-2 • • • 5o- If A: < n — 1 then there is a 0 to the left of 6 k so the identity enables the 
replacement 06k6k-i —)■ 100. If this introduces a sequence of successive ‘I’s, we apply the 
first case. Repetition of the argument, i.e. induction, completes the conversion to maximal 
normal form. 

Suppose that denotes the set of all n-digit binary sequences in maximal normal 
form for radix (j). Set Hi = {0,1}, 0.2 = {0,1, </>}. 

Recall that denotes the set of n-digit ^-integers whose leading digit is 1. For p = (j) 
one finds the disjoint union 


Vi 

= {1} 1 


V2 

= W > 

(50) 

Vn 

= (0V„_i)U(l + </.2V„_i) J 



The number of elements in is the Fibonacci number F„. 

From 0* = Un>i find 

O* = JVn 

n>l 

= Vn-1 U (1 Vn-l) 

n>l 

so 

Qn = {^Or,-l)U{l + (^^Qn-2). (51) 

The density of in Z is lim„^oo A^n+2/2” = 0. 

*** 

Let m denote Lebesgue measure on the plane induced by the Euclidean metric. Recall 
that two measurable subsets of C are essentially disjoint if the measure of their intersection 
is 0, and they are essentially identical if the measure of their intersection is equal to the 
measure of each of the sets. Thus m{A U B) = m(A) + m{B) if and only if A and B are 
essentially disjoint. 

In general, the remainder set R for a positional representation with digits A = {0,1} 
satisfies the recurrence 

pR = RU {1 + R). 

If the sets on the right side are not essentially disjoint, then they will contain a set S of 
positive measure. Now inflate 5; 5 —)■ p'^S. The set of n-digit p-integers is the subset of 
p^R whose remainders are 0. For large enough n, p'^S will contain a p-integer that, by its 
construction, will have more than one p-positional representation. This is the origin of the 
problem of non-unique representations of p-integers. 
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This description of the problem also contains the germ of a systematic alternative 
approach to finding unique representations for p-integers. 

By introducing partial remainder sets that cover R but are 
laps will have measure zero. Let us follow this approach for p 
R = [0, (j)]. Introduce Ri = [0,1] and R 2 = [0, (?!) — 1]. Then 

(f) Rl = iii U (1 + i?2) 

<p R 2 = Rl- 

Letting denote the set of n-digit ^-integers in p^Rj, these 

f^l,n = U + fl2,n-l) 

f^2,n — 1 

whence 

f^l,n = U [(f)"' ^ + ^i^n-2) , 

a relation similar to the last line of eq(|51[). In fact, both recursions produce the same set of 
i?!)-integers: eq([^ constructs integers by inserting low-order digits; eq([^ by inserting high 
order digits. That is, in any positional representation Xn-i ■ ■ - Xq for a (/)-integer, recursively 
make the substitution Oil —)• 100 so that the sequence 11 no longer occurs. 

*** 

The quadratic equation for (p implies that elements of H can be written as linear com¬ 
binations of 1 and </>. Explicitly, m<p -t- n G H if and only if m, n G U {0} and n satishes 
the condition 

(n = )V(n = 1 -|- )V((n = 2 -|- ) A (2 -|- >41 + 2 + [(m -|- l)/(/)J)), 

where [xj denotes the least integer greater than or equal to x. This condition implies that 
for a given m, n can always assume the first two of the three values { Ym/(j)\, 1 + \jn/4\ 1 2 + 
} but not any value that is not one of the three. This provides a simple way to bound 
the possible values of mcp + n. \i u) = mcj) + n € Q, then 

2 m/4 < m4 + \ rn/4\ < w < m4 + m/4 + 2 = m{24 — 1) -|- 2. (54) 

4.2 Dimension 1: Properties of radix —0 

Now let p = —4- Some integers have infinitely many representations and, more generally, 
(—(/))-integers need not have a maximal representation. Indeed, p satisfies x'^ + x = 1 which 

^"^These partial remainder sets are the ‘tiles’ for the tessellation defined by the positional representation. 


essentially disjoint]^ over- 
= 4- The remainder set is 


equations imply 


(52) 


(53) 
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is equivalent to the positional identity 110 = 001. It follows that 1 has infinitely many 
representations, viz. 

1. = no- = 11010- = ... 

Now consider the minimal representation, constructed by the substitution 110 —t- 001. 
The integer p + 1 has the representation 11 • 0 = 00 • 1, which is a remainder. It is the only 
p-integer that is also a remainder. We amend the definition to exclude polynomials p{p) 
that are equal to a remainder. Suppose p,q,r G A[X] are polynomials and the degree of q is 
m. Furthermore suppose that the radix p satisfies X^p{X) = q{X). Then the polynomial 
pip) — Qip)/p"^ is a remainder, hence not a p-integer. This is the situation we encounter 
with p = —cj). 


The (—i;f))-integers are signed. Denote the set of (—(/))-integers by X[p\ and the subset of 
n-digit (—(/>)-integers by A„[/)]. We have 

^ = E = E E (55) 

fc >0 fc >0 fc >0 

where each sum is finite. Thus A[/ 3 ] is neither bounded below nor above. 

*** 

The construction of a complete set of unique representatives of the p-integers proceeds by 
using the equivalence of digit sequences 110 = 001 to reduce representations to the least 
significant digit, i.e., to lower degree. This process is facilitated by introduction of even and 
odd p-integers. We shall say that a representation is et'en, resp. odd, if the least significant 
digit is 0, resp., 1. This notion depends on the form, but not the numerical value, of the 
representation. For instance, + p = 110 is even, but 1 = 001 is odd. 

Employing reduction to minimal normal form, set V® = {0} and, for n > 1, let V® 
denote a unique set of even n-digit p-integers with leading digit 1, and the corresponding 
set of odd p-integers. V 2 requires special attention because 11 • 0 = 00 • 1 is a remainder. 
Hence p -|- 1 ‘looks like’ a p-integer but is not one, so V 2 = 0. The following recursion is a 
constructive realization of a complete set of unique representatives for the (—()))-integers. 
For n < 3, 

V? 

whereas for n > 3, 


= 

{0}, V} = {1} \ 


= 

{p], = 0 }, 

(56) 


{p^}, Vl = {p2 + l,p2+p+l} J 


Vn 

= vOuvi \ 


n 

= (pV^i)U(p + p 2 V^ 2 ) > 

(57) 

n 

= (l + pV^i)U(l + pV%) J 
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As usual, n* = ljfc=i and Q = {0} U Q*. Set = O n (Ufc=2 ~ these are the even 
/3-integers with at most n-digits. 

The construction of is intertwined with the V^. According to the second to last line 
of eqd^, the even representations can be constructed without reference to the odd ones. 
One consequence is that this set of even p-integers can be used to define functions that 
have /9-duplication formulae. It follows that, with the initial conditions 

o? = {0}, o? = {0,/9}, nl = {0,p,p^}, 

the recursion is 

0° = (/9llti)U(/9 + /92 0t2). (58) 

The recursion summands are pairwise disjoint. Hence the number of elements in these 
sets can be calculated by observing, for example, that #(V°) = -|- #(V°_ 2 ), and 

solving the recurrence with the aid of the initial values to find #(V^) = T„_i, and then 


4.3 Dimension 2: Properties of radix 


Now consider the radix p = iy/^ with remainders A = {0,1} as above. We shall obtain 
a unique representative for each /3-integer by expressing the p-integers in terms of the (j)- 
integers, and then applying the unique expression for those integers realized by the set P 
to obtain a complete set of unique representatives for A[p]. 

The remainder set i? = {x : x = £ A, is a rectangle in C whose vertices 

are 

^ pik+l ’ ^ pik+2 ’ ^ pAk+'i ’ 

k=l ^ k=l ^ k=l ^ k=l ^ 

that is, (j) — 1, —iy/~^{2 — cj)), (2 — c/)),iy/^{2 — (j)). 

An arbitrary polynomial u) = ^kP^ of degree n — 1 in A[p] can be expressed as 


n—1 n—1 n—1 n—1 

^ ^ ^Ak+2P^^^‘^ + ^ ^Ak+d,p‘^^^^■ 

k=0 k=0 k=0 k=0 

Then 

k=l k=l k=l k=l 

= OOq + {iy/^)uJi - 4>002 - (iy/^) 

= {ujQ - (l)UJ2) + (iy/^) {uJi - (j)L 03 ). (59) 


If the ook — 4>ook+2 are restricted to a complete set of unique (—(/!))-integers, then eq(59) 
yields a unique expression for each p-integer in terms of powers of the radix p and {—4>)- 
integers. That this expression for the p-integers still requires only 2 digits for a complete 
set of positional representations relative to radix p^ is a consequence of |p^| = \ — (j)\ < 2. 
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5 Identities for partial sums of meromorphic functions 


In this section we consider partial sums of the p-Hurwitz (^-function and show that they 
satisfy certain identities even when the complete series is divergent. We begin with the 
general case, for which representations of p-integers are unique. This is the simplest, but 
also the least interesting, because it is not related to tessellations. 


5.1 If p-integer representations are unique 

Let rin = An[p] denote the set of n-digit p-integers, and assume that the representation is 
unique. Thus the coefficient of is 1. If representations of p-integers are unique, the 2^ 
n-digit polynomials in p represent distinct integers and one has the identity 

— (ptln—l) U (1 -|- plln—i). 

Introduce 

/»(»•-) = (“) 

i ’n 

A simple calculation shows that 

fn{s,pz) = fn-l{s,z) fn-l{s, Z + 1/p). (61) 


This p-duplication formula for / is an identity for finite sums so convergence is not an 
issue. 

The set of all p-integers is 

n = lim 

n^oo 

and the p-Hurwitz ("-function (p(s,z) is defined as 


(p(s,z) 


y _^_ 

o i^+y’ 


if the sum converges appropriately, whence, under the same conditions. 


= lim fn{s,z). 
n^oo 

Therefore C,p{s,z) satisfies the same p-duplication formula, viz., 

P"" Cp{s, pz) = Cp{s, z) + Cp{s, z + 1/p). 


(62) 


(63) 
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5.2 lfp = 0 

Let p = 4>. Representations of p-integers are not unique, so we must re-examine identities 
for partial sums. Introduce 


fn{s,z) 

fn{z) 



(64) 

(65) 

( 66 ) 


Since the sums are finite there is no question of convergence for any s. A calculation shows 


p''fn{s,pz) = fn-l{s,z) + ^fn- 2 {s,z/p + l/p‘^), (67) 

Pfnipz) = fn-l{z) + -fn-2iz/p+l/p‘^)--fn-2{i/p‘^)- ( 68 ) 

P P 

For each n G these are identities for finite sums. 

Now consider the limiting cases. Let Kq denote the complement of the set of poles — 
in K. In section [9TT] we show that 


fn{s,z) = lim fn{s,z) = V - -^ 

77,—^OO f 7" 4- 




{z + loY 


(69) 


exists and the series converges uniformly on compact subsets of Kq for Re (s) > 1, and that 


f{z) = lim fniz) = - + '^ 

77.—^OO r ^ ^ 




1 1 


Z + Ul UJ 


(70) 


also converges uniformly on compact subsets of Kq. As a consequence one sees that f{z) = 
cotp{z) so 


P''fis,pz) = f{s,z) + ^f{s,z/p+l/p^), (71) 

P" 

pcotp{pz) = f{z)+ -cotp{z/p + l/p^)--cotp{l/p^). (72) 

P P 

The f{s,z) are proportional to the derivatives of cotp(z). 
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5.3 If p = —0 

We shall give an exact p-duplication formula for the subset of even (—0)-integers. 
The argument follows the same lines as the previous section but is based on the dehning 
recursion eq(57) for the Put (n°)* = —{0}, = hm„_j.oo ^'^cl {^^)* = —{0}. 

Introduce the p-Hurwitz (^-function by 


and the functions 


Then 


Cp(s,z) := ^ 


fn(s,z) := ^ 


no 


1 


(Z + LU)^ 


no 


(z T 


fniz) ■— -+Y 1 




(z + UJn) UJr. 


Cp(s,z) = ^f°(s,z) 

n=l 

oo 

Cp(^) = 


n=l 


(73) 


The difference between the even (—0)-integers and the 0-integers appears from eq(58), 

Wn G {pa;n-i,p + p^a;n-2}, 


which implies 

p'fnis, pz) = /n-l(s, z) + ^fs,n- 2 iz/P + Vp) (74) 

pS 

is an exact equality, whence, if the series converges, 

p" Cp{s, pz) = Cp{s, z) + ^Cp{s, z/pZp 1/p) (75) 

is a p-duplication formula. 

This duplication formula contains an analogue of the formula ( 2 fc+i)« ~ “ 

l)C(s). Taking the limit as ^ 0, one hnds 


(/ -1) E 

(no)* 


1 

u? 


E 


1 

(1 +pw)"' 


(76) 
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6 Analytic continuation of Q(s) 


We have shown how an identity for the Riemann (^-function derived from the positional 
representation for radix p = 2 provides an analytic continuation for ('(s) to the entire 
complex plane ( eq@). Now we will follow that line of reasoning for radix p = (f>. Once 
again we emphasize that positional representations for p-integers are not unique for this 
radix. This is the interesting new fact. 

Summarizing previous notation, let denote the set of p-integers defined by the 
recursion 

Qi = { 0 , 1 } 

O 2 = {0,!,(/)} 

fin = ((?ifln-l) U (1 + (?i2fln-2) 
and put fl = limn -^00 fin- The set fin, resp. fl, is a complete set of unique representations 
for the n-digit, resp., all, p-integers. The number of elements in fin is Fn+ 2 - Let fl*, resp. 
fl*, denote the subset of non zero elements of fin, resp. fl. 

The iji-Hurwitz (^-function is 

c,i>{s,z) = ^ . 



Temporarily abbreviate Cn(s, z) = calculate 


^75.-1 


[z + Un-lY ^ {z + I + (jP'UJr,-2y 

' "77, — 


1 1 1 


^ [z!(i) p U:ri-\Y ' ^ {zj^P PijJn-2Y 

“n—1 2 

^Cn-liz/(t)+l/(t>) + ^Cn-2{s,z/(jY + l/cjY). 


Thus 


(p"" Cnis,4>z) = Cn-lis,z) + —Cn-2{s,z/(p+l/(Y), 
and the limit as n —>■ 00 yields the (/>-duplication formula for C(t>Y, z): 

(j)^ C 4 >is, 4>z) = C^{s, z) + ^CYs, zjp) + l/cY). 


(77) 


(78) 


From the defining equation for C(/>('S, 2 ), observe that hn{s, z) := 4>^ Cn(s, (pz) — Cn-i{s, z) is 
holomorphic at 2 ; = 0, so 


lim lim hn{s,z) = 


n^oo z—^co 


i)V- 


-1)4(^)> 
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whence, from eq©, 

,/.^(</-^-1)4(s) = 4(s,1/,/.2). (79) 

Now we are in a position to calculate the Laurent expansion of C(t>{s, z) at z = 0. The 
auxiliary function 


h{s,z) := C<pis,z) - ^ = Y1 


1 


■:f {z + Ljy 


is holomorphic at z = 0. Calculate 


^(•5,0) = 4(s) 


and 


ifH,, z)/dz’‘ = 


The Laurent expansion is 

4(5, 2 ;) = — + c^{s) + r(s)r(l +1) ' 


k=l 


Evaluate at z = l/(jr and use eq(79) to find 


C<t>{z) = 


i2s 


(j)2s _ _ 1 


i+E(-i) 


fc=i 


r(s + fc) C<t>{s + k) 
r(s)r(fe + i) <))2(*+fc) 


(80) 


This is the analytic continuation of C</>('2^)- We shall show below that the series converges 
absolutely and uniformly on compact subsets of the complement of the set of poles — n for 
Re (s) > 1. 

Comparison of ec 


mann ^-function, eq( 
preceding the series 


(80) with the corresponding Laurent expansion for the classical Rie- 


15) shows that they have the same general form. The factor 

lias a pole at the root of the denominator, that is, where x = (f)^ 


satisfies = x + 1, i.e., at = (j) 01 (j)^ = 1/y. Hence s = 1 is a pole. The pole of the 
corresponding factor in eq(|15[) is the root of p® = 2 with p = 2, i.e., of 2^ = 2. In both 
cases, the algebraic equation for the radix determines the pole. 


The convergence in a right half plane of the series eq(80) is enough to show that it 


provides an analytic continuation for C^l'S) to the whole plane. The series is not, however, 
useful for calculating values of the function because it converges slowly. 


7 A 1-dimensional aperiodic tessellation and the function 

cot_<^(2;) 

Now we shall turn to the connection of these functions to tessellations, and, in particular, 
to aperiodic tilings. First we focus our attention on radix p = —(p and the corresponding 
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1-dimensional tiling. The /?-integers are unbounded from below and above. Recalling the 


discussion preceding eq(58), a set Q of unique representatives of the even p-integers was 
defined by that equation. Qn satisfies the recursion 




Now dehne a cotangent function for this situation by 



(81) 


This choice of sign for ui insures that the poles of cotp( 2 ;) coincide with the p-integers. 

A duplication formula for cotp( 2 ;) can be found by repeating the argument that led to 


eq(75), or by integration of that equation with s = 1. The result is 


p cotp(pz) = cotp(z) -7 ^ cotp(z/p - 1/p^) - ^ cotp(-l/p^). 


(82) 


The p-integers are the poles, and the differences between successive poles is either cj) or 
(f)+l. The function cot^ z is approximately periodicj^on nearby inter-polar segments of the 
same length. That is, when the graph of cotp z on a vertical strip between two successive 
poles is shifted onto a congruent such strip, the graphs of the function are approximately 
equal. Moreover, the graphs on the strip of width cj) are approximately equal to the graphs 
on the strips of width (p + 1 after they are simultaneously scaled by a fixed factor. The 
precise relationship is given by the p-duplication equation. 

If the poles of cotp(z) are interpreted as the endpoints of the intervals they delimit, these 
intervals are the tiles of an aperiodic inflationary 1-dimensional tessellation. Hence there 
are two types of tiles: intervals of length cj), and intervals of length 4> + 1. The poles, and 
therefore these tiles, are distributed aperiodically. Figure ^ displays the graph of cotp(x) 
and part of the aperiodic 1-dimensional tessellation formed by the even p-integers. 


^®The phrase “approximately periodic” is used in an informal sense. 
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Figure 3: Graph of cotp(x) for the subset of even integers for radix p = —cj), showing the 
location (indicated by the red bars) of the poles at the even p-integers. The successive 
intervals constitute an aperiodic tessellation of K. 
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8 Function theory on inflationary tessellations. The special 
tessellation 

8.1 Preparations 

Now we turn our attention to function theory on inflationary tessellations of the plane. 

Consider an aperiodic inflationary tessellation of C constructed from n types of tiles 
defined with respect to the equivalence relation of congruence. Since copies of a tile of a 
given type are congruent, there is a 1-1 correspondence of their points 

The special tessellation offers an interesting example. It is ‘special’ because it is aperi¬ 
odic; it consists of four tile types, all of which are rectangles; three of the four tile types are 
similar. The radix (the ‘multiplier’) is p = iy/~^ where (j) = is the golden number; the 
tiles in the tessellation are rotated by the multiple of 7r/2 that is the power of p is associated 
with in the positional representation, but they are never reflected. This tiling provides a 
convenient test laboratory for studying functions defined on an aperiodic tessellation. 

Consistent with previous notation, A = {0,1} is the set of digits, and A[/5] denotes the 
set of polynomials in p with coefficients in A; these are the ‘integers’ relative to radix p. 
There are repetitions. We shall construct a set of unique representations for a subset of 
the p^-integers. Like the set of ‘even’ (—()))-integers considered in section]^ is a set of 
unique representations, but it is not complete. Then we will use 11 to construct a set of 
unique representations for a subset of the p-integers. 

p is the largest imaginary root of = 1. As we know (because p^ = —<(>), this 

implies that infinitely many elements of A[p] have more than one expression in A[p]. For 
instance, 

m 

1 = p4’« + ^p4fc-2, rn>l (83) 

k=l 

Division by p^ shows that a sum of p-integers can be a remainder: 1/p^ = 1 -|- p^; indeed, 
p-2 = -l/<^ = 1 — (j) — 0.61 has infinitely many representations as a sum of non-negative 

powers of p. The various realizations of 1 -|- p^ have the least absolute value of any element 
of A[p]. 

*** 

Now turn to the associated tessellation. The remainder set is 

i? = jz : z = ^2;fcP“^| , Zk€A, (84) 

that is, the rectangle 

R= {x + iy : p^/(j)<x < !/(/>, p^/(j)<y < p/cj)] . 

^®Up to a set of measure 0. 
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An arbitrary point of C will be covered by a sequence of inflations of R because 0 is an 
interior point of R so the n-th inflation contains the disk of radius r” centered 0 if r < !/(/) 
because the disk of radius r centered at 0 is contained in R. It follows that an essentially 
disjoint decomposition of R that behaves well with respect to inflation by p will lead to a 
tessellation of the plane. We will construct such a decomposition by the method introduced 
in the previous section. 

Although R contains the origin as an interior point, its successive inflations cover the 
plane but do not tile it because the decomposition R = (1 + R) U R is not essentially 
disjoint. The translated copies overlap each other. Following the procedure previously used 
in 1-dimension |3]. R can be written as the essentially disjoint union of 4 subsets Ri that 
satisfy eq@; cp. fig. 0 The subsets are the parts of R contained in the quadrants of the 
plane. This situation is described by the tessellation dehned by 

p Ri = i?2 

pR2 = Rs 

pRs = (1 + R 3 ) U Ri 

p Ri = (1-l-i?2) U 



Figure 4: Radix p = Decomposition of the remainder set R into an essentially disjoint 
union of remainder sets Ri, R2, R3, Ri- Ri is the (magenta colored) rectangle in the first 
quadrant. 



The tile types Ri, R 2 and iia are similar. i?i is the intersection of the remainder set R 
with the first quadrant. This will be called the special tessellation. 

A portion of the 10-digit decomposition of Ri is shown in fig.[^ Note that complex 
conjugation does not appear in the equations which implies that the tiles in the decompo¬ 
sition equations are not reflected. The special tessellation is aperiodic. 
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Figure 5: Radix p = i^/^. 10-digit tessellation of remainder set Ri. See text for details. 



Proof: The proof by irrationalitjj^ applies. The partition matrix, 

/ 0 1 0 0 \ 

0 0 10 

^ ~ 0 0 11 

V 1 1 0 0 / 

has largest real eigenvalue f so p = yf—fj) is a multiplier (Note that U is the 

partition matrix for a 2-dimensional tiling, so the eigenvalue is the square of the absolute 
value of the inflation factor. The tile type densities, which can be calculated from limits of 
quotients of the entries of t/”, are irrational, so the tiling cannot be periodic; cp. [1]. □ 


8.2 Pp{z) and its p-duplication formula 


According to eq(59), each p-integer can be written in terms of (p^)-integers, that is, (—</>)- 
integers, in the form 

UJ = uj' + puj". 


Suppose that Re (s) is sufficiently large. Let 11 denote a set of unique representations 
for the p^-integers, and Vln C 11 the representations of the n-digit p^-integers. Note that 
p^ = —cj). The particular set H of unique p^-integers, i.e., (—(/))-integers, that we have in 
mind is dehned by the following relations: 


Hi = {0,1} 

II 2 = {0, !,-</)} 

Hn = (-0) U (1 -F (-0)^ 11„_2) ^ 

i"cp. a. 


(86) 
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Then O = U^i = linin->.oo ^n- Put — {0}. Induction shows that #0,n = Fn+ 2 - 

Define the p-Weierstrass p-function Pp{z) and the related p-Hurwitz function p^p\z) 


b} 


18 


Pp{z) = 

a;eSl’ 

7 - Re(s)>2. 

^ {z-ujY' ^ ^ 


\{z — 0jY {—^Y 
u)&n* ^ ^ ^ ' 


4''(^) = 




Pp”^(z) is proportional to the n-th derivative of pp{z). 
In terms of D, the p-Weierstrafi function is 


“ 72 + X] 


I ^ f 1 


^ \z — {uj'+ p ijj"Y [uj' + puj"Y 


with a corresponding expression for p^Y\Y- We shall show that these sums converge in 
section 19.31 

Figures and [^display the modulus of pp{z) and p'p{Y respectively. 


(87) 

( 88 ) 


(89) 




^®Notice that the sign convention departs from what is customary for the Hurwitz (^-function, in order 
to position the poles of the meromorphic function at the p-integer points of the associated tessellation. 
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Figure 6 ; Radix p = The modulus of the p-Weierstrafi p-function. 6 digit series 

expansion of the function. The domain is the remainder set R. The decomposition of R 
into four essentially disjoint tiles i?i U R 2 U R 3 U R 4 is evident (Cp. figlU page[^. 


Figure 7: Radix p = The modulus of dpp{z)/dz. 6 digit series expansion of the 

function. The domain is the remainder set R. The decomposition of R into four essentially 
disjoint tiles Ri U i ?2 U R 3 U R 4 is evident (Cp. figj^ page 35). 
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In constructing the p-Weierstrafi function and deriving its ^-duplication formula we 

shall make use of helping functions Pp,ni,n 2 i^) Pp%i,n 2 {z) that are finite sums whose 

As), 


limits will be pp (z). Introduce the finite sum 

(s) ( . ^_ 

yp,ni,n2\'^J ■ Z_^ / / , I 'iV 

l^ni +P‘^n2)j 


(90) 




and the corresponding finite sum for Pp,ni,n 2 {z)- If the limit function exists, we have 

712^00 

with the corresponding limit formula for Pp,ni,n 2 {z)- 

The following identity for finite sums is the p-duplication formula for pp%,n(z): 

1 (s) f 1 (s) f Z 1 \ 

P Pp,n,niP^) — —Pp^n-l,n-l + ^Pp,n-l,n-2 ~ ^ 


(91) 


n 3 s 


o'^p,n—2,71—2 \ 3 


i + p 


1 




1 (s) 

4s Pp,n—2,n—3 \ 4 


P" P" J P- 

Z ^ + P 


-,3s Pp,n—2,n—2 1 


(92) 


Proof: The proof consist of a verification that the various terms on the left side of eq(92) 
can be re-written as the terms on the right side. This is accomplished by using eq(86). The 
general summand on the left is 


/ a;; + p< 

V P 


—s 


From eq(86) we see that one way Un can be expressed is ujn-i, Apply this to to 

find 

T P II , n \ 

--^ P i^n-l + P‘^n-l) ) 


and move the factor p outside the sum to obtain p ‘^p^p\_in_i {z/p), which is the first 
term in eq(|92[). The other terms are handled similarly. □ 
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If the series converges, the p-duplication formula for (z) is 


- 1I ±_±\< I ± _ 


P'ppipz) = —P^p 


+ 


1 


,2s 


+ 


1 


,3s ^P 




o3s ^P 


,4s ^P 


1 + p 


+ 


(93) 


Proof: Take the limit as n —>■ oo in the previous theorem. It will converge, as we shall 

later see, for sufficiently large Re (s). □ 


The derivation of eq(8.2) did not depend on s being an integer; it sufficed that Re(s) 
be large enough for the underlying series to converge. Therefore, the p-duplication formula 
is valid for the p-Hurwitz ("-function for arbitrary s G C with large enough real part. We 
shall use this fact below. 


Regarding pp{z) itself, application of the same method yields: 


The p-duplieation formula for pp{z) is 



9 Convergence properties 

It is time to take up the question of convergence of the infinite series we have been freely 
manipulating. In the framework of this paper, we shall not be concerned with establishing 
the tightest conditions for convergence; rather, we merely want to show that the series 
converge for some potentially significant values of the variables so that the functions we 
have introduced, and the identities they satisfy, make sense. 

The convergence problems concern the p-Riemann (-function; the p-Hurwitz (-function; 
and the p-Weierstrafi p-function. In the classical case one first proves convergence of the 
series that defines the Riemann (-function, reduces the convergence of the Hurwitz (- 
function to it; and applies the same approach to the Weierstrafi p-function, which can be 
thought of as a 2-dimensional version of the Hurwitz function. 
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The new series are sums over a set of p-integers. The convergence proofs follow the 
general pattern of the corresponding classical case; cp., for instance, the famous textbook 
of Whittaker and Watson |8] . The principal difference is that the summands are collected 
according to the leading term in the positional representation of the p-integer summation 
variable. 

As above, let 1 < |p| < 2 be a radix and A = {0,1}, and let A[p] denote the set of p- 
integers. In general, the notation Yls* stands for the sum over the index set S* := S — {0}. 
The convergence of the series will be subdivided into cases, according to the nature of p: 
p G K'*'; p G and p G C — Within each case we shall also have to consider whether 
each p-integer has a unique positional representation. 

The idea behind the proof of convergence is to collect terms in the sum that have the 
same highest power of p. Let A[p]„ denote the set of n-digit p-integers. The set of p-integers 
whose highest power of the radix is p"' is Vn := A[p]„ — A[p]„_i. 

Begin by considering the p-Riemann (^-function; cp. eq([n]). 


Cp(s) := Yj 

a;eA[p]* 


1 


9.1 p G 


This is the simplest case because a p-integer is a sum of non-negative monomials. 

Let ujn = p””^ -|- Ylk=o ^kP^ G Vn be a p-integer whose highest power of the radix is n. 
Then 




n=l V„ 


1 




and 


Evidently 


oo .. oo .. 


n=l V„ 


p" ^ < Wn < 


n=l V„ 

p" - 1 
p- 1 ■ 


(^)- 


Since the 5^ are 0 or 1, there are at most 2*^ ^ distinct n-digit p-integers, whence 


P-^ 


< 2^-^— < 


1 


in—1 


p" - 1 


Wn P"- ^ ' 


SO 


in—1 




oo ^ 2 ^ 


p" 


1 


1 - 2p-^’ 


(96) 


n=l^ n=l 

so we have proved that if 1 < p < 2 and A = {0,1}, then Cp{s) converges absolutely for 
Re(s) > log 2/log p. 
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Note that Re (s) > 1 for all radices. Recall that the representation of p-integers need 
not be nnique. Nevertheless, the lemma applies, for instance, to p = (p where positional 
representation of p-integers is not nniqne because the radix is positive. In this case, we 
have convergence for Re (s) > log 2/ log (/> ~= 1.44 .... 


The radix p = 4> generates more than one representation for some p-integers. As we show 
in section 4.1, the cardinality of Hn, i-e. the number of pairwise unequal representations 
for n-digit (/>-integers is Fn +2 where Fn is the Fibonacci number. The convergence criterion 
for 


Cpis) = Y, 


U! 




is derived as above, but now the cardinality of is Fn rather than 2 "' so |CPC'S)! < 
0 (n-pRe (s) ■ Since Fn is the closest integer to the series converges when 


^ n 

fARe{s)\n-l 
n>l ' 





n—1 


converges, that is, when Re (s) > 1. This improves on the general condition Re (s) > 
log 2 / log p given above. 




Convergence of the p-Hurwitz is dealt with in the usual way, by noting that 


E 


1 

{z + uY 


< 


E 


1 


Z -huj 


Re (s) ’ 


and that for z in a compact set, |2 + a;|, apart from finitely many terms, the terms of 
the series are dominated by the corresponding terms of Cp(s)) foi' which convergence has 
already been shown. 


9.2 p e K“ 

Now suppose that — 2 < p < — 1. The p-integers have negative as well as positive values 
that are unbounded as the number of digits is increased. Again consider an n-digit p-integer 
whose leading coefficient is 1 . 

Wn = ^2kP^^ + P Yj ^2fc+lP^^, 

fc=0 k=Q 
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where the sums are non negative. Un is positive, resp., negative, if n 
If n — 1 = 2Z, then 


p‘^^ + P 



— ^ 


p2l+2 _ 1 
- 1 


1 is even, resp., odd. 


A p-integer that is closest to the origin is p + 1. If p / —2, it is unique. Let ojn denote, as 
above, a p-integer that is expressed by exactly n-digits, so that the coefficient of p”“^ is 1. 

For p < — 1 we cannot hope to bound \ujn\ away from 0 by a bound that increases with 
n in general. For example, if p = — ( ), then 


k=l 

SO the set of p-integers with an odd number of digits always contains 1. This implies that 
sums such as X]a[p]* have an unlimited number of terms equal to 1“^ and cannot 

converge. This shows that a radix for which the representation of a p-integer is not unique 
require special, and perhaps individual, consideration to eliminate duplications. 

We restrict consideration to p = —cp and to summations over the set of even p- 
integers defined by eq(57). Recall that #12° = Each a; G 12° has the form 


n—1 


U = 


Y.^kP^- 


k=0 


Denote the subset of elements of 12° whose leading coefficient 5n-i = 1 by V°. Thus 
V° = 12° — f2°_^. Evidently #V° = F„-i. The non-zero elements of V° are negative, resp., 
positive, if n is even, resp., odd. 

We seek a lower bound for Un G V°, and find < \oJn\- Thus 


V - 


< 


E 


1 


0 ;“ 


EE 

n=l (VO)* 
oo 

EE 

rx=l(VO)* 

OO 

E 


0 ;“ 


(rn-i)‘ 




\S-l 


E 

71=1 


1 


43 











whence the series converges for Re (s) > 1 + 




3. 


9.3 p G C — K with A = {0,1} 


We confine ourselves to the complex radices that are the classical complex multiplications 
for elliptic functions, and to p = 

In the classical cases, where p is an imaginary quadratic integer, an arbitrary p integer 
has the form m+np, where m,n G z|^For these cases, the problem of convergence proceeds 
as usual: a sequence of artfully arranged similar nested parallelograms covers the non-zero 
p-integers; a counting argument bounds the number terms on each parallelogram; and an 
estimate for a lower bound is found that is proportional to the index of the parallelogram. 

In the general 2-dimensional case we have neither the 2-dimensional basis nor the 
parallelograms. Instead, we partition 12* into sets as we have on other occasions. In order 
to do this we must make explicit the structure of a complete set of unique representatives 
n for the p-integers. 

For p = the geometry is much simpler, and this is why we use it as an example. 
The p-integers are arranged on a sequence of concentric rectangles centered on the origin. 
They are not distributed in the regular manner found in the classical case, but they are 
nevertheless easy to count. 

Number successive similar rectangles starting with 1, and denote the A:-th rectangle by 
rectfc. The number of p-integers on rect is 8k (cp. hgj^ and the total number of p-integers, 
including 0 , in the closure of the fc-th rectangle is 1 - 1-8 Yln=i ^ ~ (2A:-|-1)^. These rectangles 
are successive inflations of a remainder set for radix p, which provides the foundation for 
a detailed proof. 

The minimum distance from the origin to rect^ is achieved along the positive real axis. 
The differences between successive real p-integers is either 1 or cf), whence the distance of 
points on rectfc from 0 is greater than or equal to k. This will suffice to show that the 
various series converge for interesting values of the variable. 

Consider, for instance, the p-Riemann (^-function, 




We find 


\Cpis)\ = 


V- 


oo 


oo 

SEE 

k=l rectfe 


\^k I 


1 8k 

Re (s) / ^ j^Ke (s) 

k=l 


E 

k=l 


1 


/jRe(s)-l ■ 


^®The p-Riemann (^-function (p(s) = cj can be interpreted as a special value of an Eisenstein series 
when s G 2Z~^. 
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which converges for Re (s) > 2. This simple result is enough to conclude that the series for 
the derivatives of the p-Weierstrafi pp-function converge, because they are just the numbers 
Cp{k + 2), k G Z~*~. 

Regarding pp{z) itself, we follow the standard proof. Suppose that z is an element of a 
compact subset of C — 12. The absolute value of the general summand satisfies 

1 1 

(z —w)2 OJ^ 

whence the series for pp{z) is uniformly convergent on compact subsets of C — 12. 


\z — 2ijj\ 

'I 191 19 

\z — up\ojp 

1 


9.4 Analytic continuations 

Equation (15) and its special case eq@ for p = 2, and eq([80l) for p = (p, provide an analytic 
continuation of the respective Cp(s) to C. Comparison shows that the two expressions have a 
similar form, apart from a common factor that depends only on the radix and s, the power 
of p that appears in the denominator of each summand, and the intrinsic difference between 
the two p-Riemann (^-functions. Here we shall investigate convergence of the shared series 


k>l 


r(s + kp Cp('® T 
r(s)r(A: + l) pds+fc) ’ 


(97) 


where c = 1 or c = 2. 

We have shown that C,p{s) converges for Re(s) sufficiently large, and it is evident that 
for any e > 0, there is an such that |Cp(s) — 1| < e; hence, as far as convergence is 
concerned, we can ignore the factor ^p(s + k). 

Temporarily write s for Re (s) to simplify notation. Application of Stirling’s formula 
shows that 

r{s + k) 1 (l + s/fc)^+i k^-^ 

r(s)r(A; + l) |p|d^+fc) "" ^ (1 + 1/A:)2 \p\cis+k) ~ |p|c(s+fc)- 

The terms go to 0 with increasing k. The comparison test shows that the ratio of successive 
terms is ultimately, i.e., for sufficiently large Re(s), l/|p|^ < l/|p| < 1, so the series 
converges to the right of the pole give by the factor preceding the series. This factor was 
/{(p^ — 1) if the representative polynomials A[p] are unique, and — 1) for 

p = p. 
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10 Remarks 


10.1 About the color coding of the graph of a complex valued function 

The graph of the modulus f{z) of a complex valued function is colored to show the value 
of the argument of f{z). The colors are drawn from the circle of hues so that positive real 
numbers (argument = 0) correspond to red. As the argument increases from 0 to 2tt, the 
hues range from red through yellow, green, teal (negative reals) blue, and magenta. 


10.2 About the numerical calculation of identities 


I have checked many of the new - especially the duplication formulae - numerically. The 
calculations are numerically unstable because they require cancellation of values in neigh¬ 
borhoods of the many poles of the functions. A few remarks about how this can be done 
may be helpful. 

The problem can be illustrated by examining expressions like f{s,z) = (z + uj)~^ that 
appear in many of the formulae. The p-duplication formulae involve expressions such as 

P‘'f{s,pz), f{s,z), -f{s,z/p+l/p^). 

P 


which the theorem tells us cancel. 

Because all numerical calculations are performed with a fixed maximum number of 
digits, it is essential to insure that only terms with the same number of digits appear. Recall 
that Qn denotes the set of n-digit p-integers for some radix p. Set fn{s,z) = (z + ujn)~^ 
where ujn G fin • Multiplication and division of z by powers of the radix changes the number 
of digits. For instance. 


P"" fn{s,pz) 


1 

{z + OJn/pY 


SO the numbers iVn/p for which the low order digit of cOn is zero actually belong to fln-i 
rather than fin • Identities such as those in section insure that the terms in a calculation 
are effectively of the same order, thereby enabling cancellations to an arbitrarily high degree 
of accuracy. 
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